Scaling Limit of Vicious Walkers, Schur Function, 
and Gaussian Random Matrix Ensemble 
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We consider the diffusion scaling limit of the vicious walkers and derive the time-dependent spatial- 
distribution function of walkers. The dependence on initial configurations of walkers is generally 
described by using the symmetric polynomials called the Schur functions. In the special case in the 
scaling limit that all walkers are started from the origin, the probability density is simplified and it 
shows that the positions of walkers on the real axis at time one is identically distributed with the 
eigenvalues of random matrices in the Gaussian orthogonal ensemble. Since the diffusion scaling 
limit makes the vicious walkers converge to the nonintersecting Brownian motions in distribution, 
the present study will provide a new method to analyze intersection problems of Brownian motions 
in one-dimension. 
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The problem of vicious walkers was introduced by 
Fisher and its application to various wetting and melt- 
ing phenomena were described in his Boltzmann medal 
lecture Recently, by using the standard one-to- 

one correspondence between walks and Young tableaux, 
Guttmann et al. Q and Krattenthaler et al. showed 
that the exact solutions for some enumeration problems 
of vicious walks |4 ^ are derived from the theory of sym- 
metric functions |B[ or the representation theory of clas- 
sical groups . Important analogies between the ensem- 
bles of Young tableaux and those of Gaussian random 
matrices were reported by Johansson and Baik ||] 
and Nagao and Forrester jl^ studied the vicious walker 
problem using the random matrix theory of the Gaussian 
orthogonal ensemble (GOE) [l2|. 

The purpose of this Letter is to demonstrate more ex- 
plicit relations among the vicious walker problem, the 
Schur function [|l3j, and the GOE, by performing the 
diffusion scaling limit of the vicious walkers. We derive 
the time-dependent spatial-distribution function of walk- 
ers in this scaling limit. The dependence on the initial 
configurations is generally described by using the Schur 
functions. We show that the case, in which all walkers 
are started from the origin, can be treated, and in this 
special case the probability density of positions of walk- 
ers at time t = 1 is identified with the probability density 
of eigenvalues in the GOE. Since the distribution of ran- 
dom walkers converges to that of Brownian motions in 
the diffusion scaling limit, the present analysis will solve 
some intersection problems of one-dimensional Brownian 
motions. More applications to the probability theory will 
be reported elsewhere 

Vicious walks are defined as a subset of the sim- 
ple random walks as follows. Let {-Rfe'}fe>o,* Sz In = 



{1,2,- • ■ ,n}, be the n independent symmetric simple ran- 
dom walks on Z = {• • • , —2, —1, 0, 1, 2 • • •} started from 
n distinct positions, 2si < 2s2 < ■ ■ • < 2s„, Si € Z. Fix 
the time interval m as a positive even number. The total 
number of walks is 2™", all of which are assumed to be 
realized with equal probability 2~™". Now we consider 
a subset of walks such that any of walkers does not meet 
other walkers up to time m, that is, the condition 



Rl} <Rt' <■■■< Rl 



1 < Vfc < m 



(1) 



is imposed. Such a subset of walks is called the vi- 
cious walks (up to time m) Let Nn{m; {ei}\{si}) 
be the total number of the vicious walks, in which the n 
walkers arrive at the positions 2ei < 2e2 < • ■ • < 2e„ 
at time m. Then the probability that such vicious 
walks with fixed end-points are realized in all possi- 
ble random walks started from the given initial config- 
uration is iV„(m; {ei}|{si})/2™", which is denoted by 
Vn{{Rk}T=o''Rm = 2ei) in this Letter. We also define 



K({i?^'}r=o)- 



E 

ei<e2<-<e„ 



Recently Krattenthaler et al. ^ evaluated the asymp- 
totes for large m of ^({-Rfe'jfcLo) the two special 
initial-configurations, (i) Si — i — 1 and (ii) Si — 2{i ~ 1), 
as 

K({i?fe'}r=o) = (1 + 0(l/m)) , 

where 

n/2 

(2"/7r)"/4 Yli2i - 2)! if n = even 

(2"+V7r)(""i)/4 Yl (2j-l)! ifn=-odd, 
1=1 
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with 6„({i - 1}) = 1, 6„({2(i - 1)}) = 2"("-i)/2. 

We found that their result can be immediately gener- 
alized as 



6„({s(i - 1)}) = s"(«-i)/2 for any s = 1, 2, ■ 



(3) 



This observation suggests that we can take the scaling 
limit such that L oo with the time interval m — Lt 
and the initial spacing of walkers s = where t is 

finite. 

In this Letter we consider the diffusion scaling limit] 
setting TO — Lt,Si — ^/Lxi/2, Ci = \/Lyi/2,i G /„, and 
taking the limit i — > cxd. The key lemma, which will 
be proved shortly, is the following. For given t > and 
< xi < X2 < ■ ■ ■ < Xn (xi € Z),yi < y2 < ■ ■ ■ < yn, let 
^(x) = {£,i{x), ■ ■ ■ ,^ri{x)) be a partition specified by the 
starting positions {xi} as 



then 



^i{x) = Xn-i+i - (n-i), i e In, 




(4) 



X exp 



), (5) 



l<i<j<n 



where • • • , z„) is the Schur function jl^. We 

consider the rescaled one-dimensional lattice Z/(VX/2), 
where the unit length is and let denote 

the symmetric simple random walk starting from x on 
Z/(Vi/2). Then (|) implies that 



lim Vn ( f i?fc' 1^* ; RVt € [y., y^ + dy^ 
-/n(t;{yj|{xj)d"2/, 

where {yi}|{a;i}) is defined by the RHS of (||). 
In order to normalize /„, we consider the integral 



(27ri)"/2n! 



d''y fn{f{y^}\{x,}) 

yi<-<Vn 



xe-E«.V2t -Q |eK/t_eWt|^ (6) 

l<2<j<n 

where we have used the fact that with the absolute values 
the integrand is invariant under permutation of yt , since 
the Schur function is a symmetric function. The spatial- 
distribution function in the scaling limit of the n vicious 
walkers at finite time t is then given by dfin — fnd^y/Mm 
or more explicitly 

d^ln{t■, {yi}\{xi\) = gn(t; {yi}\{x.,})d"'y, 



with the probability density 

i{yi < y2 < ■ ■ ■ < yn} 



9n(t; {y^}\{x^}) = 



Zn{t; {Xi}) 



xs^(,)(e^i/*,...,e^"/*) 



X exp 




1<2<J <n 



) (7) 



for < a;i < • • • < Xn,Xi G Z, where is 1 if 

uj is satisfied and is zero otherwise, and Zn{t;{xi}) = 

(27ri)"/2eS'''/^*Mi(*;{a:»}). It should be noted that, if 
the initial configuration is a;^ — i — l,i G /„ , then 
£,i{x) = and the Schur function in (||) and is 



= 1. 



Now we give a proof of 
square lattice Z^, 



Define a subset of the 



= {(a;, y) eZ'^ : x + y ^ even, < y < to}, 

and £rn be the set of all edges which connect the nearest- 
neighbor pairs of vertices in Cm ■ Then each walk of the 
i-th walker, i € In, can be represented as a sequence 
of successive edges connecting vertices Si = (2si,0) and 
Ei — (2ei,TO) on [Cm,£m), which we call the lattice path 
running from Si to Ei. If such lattice paths share a com- 
mon vertex, they are said to intersect. Under the vicious 
walk condition (]l|), what we consider is a set of all n- 
tuples of nonintersecting paths [p^ . Let Tr{S — > E) be the 
set of all lattice paths from S to E, and 7ro({S'i}f^]^ 
{Ei}f^i) be the set of all rt-tuples (tti, • • • , 7r„) of nonin- 
tersecting lattice paths, in which tt; runs from Si to Ei, 
i € In- If we write the number of elements in a set A 
as \Al then 7V„(to; {ej|{s,}) = |7ro({5JJLi ^ i^^Kk)! 
and the Lindstrom-Gessel-Viennot theorem gives p^ . 



Nnim;{e,}\{si}) 



, det {\n{S, 



Since \tt{Sj E,)\ 
determinant 



we have the binomial 



= 2-""' det 

l<i,j<n 



TO 

m/2 + Sj 



Application of Stirling's formula yields 



lim 2-'^*"(\/I/2)" det I I , ^, 

L^oo l<i,j<n\\Lt/2+ VL{xj 



Lt 



det hm 2 

1<2,J<71 \L — >00 



Lt 



^^^\Lt/2 + VL{x, 
(27rt)-"/2g-^(x?+2/n/2t det fe^^-^' 

Ki.j <7i \ 



- y^)/2. 
yO/2 



= det 

l<i,j <n 
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We rewrite the determinant as 



det f e^^ 

l<i,j<n \ J 



^ deti<,j<„ ((e^./*)^^-.+i 
deti<,j<„ ((e!'./*)"-J) 

where A„({zi}) is the Vandermonde determinant 
A„({z.})^ det (zr')= n (^^--^O- 

l<i<j<ri 

Combining with (Q) and the definition of Schur function 
[p^ completes the proof of (||). 
Since 

= S4(a;)(e^i/*,---,e«"/*)A({e^'}) = det fe^^'^"/*) 

l<ij<n V / 

for any integer where Ix ~ {ixi, • • • , we can prove 
the foUowing scaling property for the scahng-hmit prob- 
abihty density. For any integer i 

g^it- {yi\\{xi\) = rgnilh- (8) 

Using this property, we can generahze the expression (j^) 
for any rational numbers, xi < ■ ■ ■ < Xm and then using 
the connectedness of real numbers, for any real numbers 

{Xi}. 

Next we study the t ^ oo asymptotes of the above 
results. Since Q 

hm s^(,)(e^i/*,...,e^"/*) = s^(,)(l,l,...,l) 

= n ^^^u{x^}) 



l<'i<j<n 



J ~ « 



(remark that this definition of 6„({x}) is consistent with 
(|), and 

the normalization factor is asymptotically 
Mnit;{x^})^t-""/^il + Oil/t)) 

{2ny 

= i-"("-i)/4(l + C'(l/t)) 

bni{Xt}) 



l<i<j<n 



{2n) 



l<«<j<n 



where Ui = yi/\/t. The last integral is the special case 
(7 = 1/2 and a = 1/2) of 



J] 
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l<'i<_;<n 



"/2('9«^-"('y(n-l) + l)/2 fr r(l + »7) 



(27r)"/2(2a) 



r(i + 7)' 



which was derived by Mehta (eq. (17.6.7) on page 354 in 
JlH ) as a consequence of Selberg's integral. Here T{x) is 
the Gamma function with the values r(3/2) = \fT^ 1'^ and 
nr=i r(l + V2) = 2-"("-i)/2(V^/2)"n! a„, where a„ is 
given by (^. Then we have 

= t-"("-i)/^2-"("-i)/2a„6„({xJ) (1 + 0(l/<)) (9) 

as i tends to infinity. Using (^ the asymptotic form in 
t — + cxD of the probability density is given as 

<?n(i; {V^}\{X^}) = l{2/l < 2/2 < ' ' ' < 2/«}c„t-"("+l)/4 

xe-i:^'.V2t -Q (l + 0(l/t)) 

l<i<j'<?T, 

with c„ = 2"("-i)/V{(27r)"/2a„} and thus 

lim d^„(t; {\/tjTUi}\{xi}) = < U2 < • • • < u„} 
(c„r-"("+i)/4 e-E"?/2- J] (u, -u,)d"ii (10) 

l<i<j<n 



t — >oo 



for any r > 0. We note that the diffusion scaling limit, 
L — !■ oo with m = Lt,Si = y/Lxi/2^ei = \/Lyi/2, and 
the above t ^ oo limit with yi = yjt/rui are combined 
to define a limit L' = Ltjr ^ oo with m = L't, Cj = 
\flJuil2 and Si in the smaller order than \JU . Since 
Sj/ei — > in this hmit, the limit ( |l0|) should be regarded 
as the scaling limit d/i„(t; {yi}|{a;i}) for xi = 0, i S 
That is, when all walkers are started from the origin, the 
probability density in the scaling limit should be 

5«(t; {yj|{0}) = i{yi < y2 < • • • < yn} c„i-"("+i)/4 



X exp 



1 " \ 



(11) 



l<i<j<n 



By using the scaling property (|8|), we can give another 
explanation for the reason why the above limit procedure 
gives the case for Xi = as follows. The scaling property 
(^ can be written for distribution functions as 

d^in{t■Av^}\{x^}) = d/i„ (^2^; {%} | {£x,}) 

for any integer I. Then, if we set Xi = Q and T — £'^t, we 
have the invariance 

d^in{t■, {y^}m) = d^i,,{T;{^t y,}\{0}) (12) 

for any T > 0. Then we can take the limit T —^ oo, which 
will derive the expression (|Tl|). 

The result (111]) is very important, since if we set t = 1 
and assume that yi < ■ ■ ■ < y-n, then we have the equality 

g„(l;M|{0}) = n! ^^^^({yj) 

where 3,^°^({2/i}) is the probability density of eigenval- 
ues of random matrices in the GOE pTl]. The equality 
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( p^ with T = 1 implies that {yi}|{0}) is a func- 

tion of n variables {ui = jji/Vt}- In other words, the 
dynamical scaling with the classical exponent 6 = 2 (see, 
for example, chapter 16 in is established, in which 
the scaling function is exactly and explicitly given by 
nl X with the scaling variables Ui = iji/t^^^. 

In the present diffusion scaling limit, the random walks 
converge to the Brownian motions in distribution. Then 
gives the asymptote in t ^ cx) of the probability 
that the n independent one-dimensional Brownian mo- 
tions started from xi < X2 < ■ ■ ■ < Xn {xi & 'Z) do 
not intersect up to time t. The scaling property (||) is 
nothing but the diffusion scaling of Brownian motions, 
J^ni^^t', {^Xi}) = Afn{t;{xi}). Moreover, we can define 
a diffusion process by giving the transition probability 
density as 



Pn{s; {yi}\{x^}) = lim 



fn{s;{y^}\{x^})^^n{t;{y^}) 



Afn{t + S; {X^}) 



b„{{xi}) 



fn{s;{yi}\{x^}), 



where (^) has been used |T^. Isozaki and Yoshida 
studied the n = 2 case as the limit process of the two 
friendly walkers in the "wetting phase" . 

In summary, we have performed the diffusion scaling 
limit of the vicious walkers and determined the time- 
dependent spatial-distribution functions in the limit. 
There the Schur function plays an important role to de- 
scribe the initial-configuration dependence. The relation 
with the Gaussian orthogonal ensemble of random ma- 
trices was clarified. In the present study we have con- 
structed new continuous-time processes on the real axis 
for arbitrary numbers of walkers n. We remark that the 
results reported in this Letter may hold also in the fully 
"wetting phases" of the models of friendly walkers 
pO| , which are related with a famous un- 

solved problem, the directed percolation problem. Fur- 
ther study of the distribution functions and diffusion pro- 
cesses derived in this Letter will be desired. 
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